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Coherence  theory  is  used  to  analyze  the  statistical  properties  of  ocean-acoustic  intensity  fluctuations 
measured  after  saturated  multipath  propagation.  Previous  analyses  in  this  area  have  been  implicitly 
limited  to  certain  special  cases  for  which  the  time-bandwidth  product  of  the  field  received  from  a 
given  source  is  unity.  In  this  paper,  the  statistical  description  is  extended  and  generalized  to  be  a 
function  of  measurement  time  and  temporal  coherence.  As  a  result,  the  well  known  5.6-dB 
transmission  loss  (TL)  standard  deviation  of  Dyer  is  found  to  be  a  special  case  of  a  more  general  TL 
standard  deviation  that  approximates  4.34  sqrt(l//r)  dB  when  the  time -bandwidth  product  /jl  is  large. 
Therefore,  the  TL  standard  deviation  approaches  zero  for  increasing  j±,  as  it  must  in  the 
deterministic  limit  of  an  arbitrarily  large  sample  size.  A  similar  generalization  is  obtained  for  the  TL 
mean,  from  which  it  is  found  that  the  sonar  equation  must  be  corrected  for  a  //-dependent  bias  that 
vanishes  in  the  deterministic  limit  of  large  //.  Additionally,  asymptotic  analysis  shows  that  intensity 
statistics  in  the  saturated  region  converge  to  a  log-normal  distribution,  where  fi> 4  is  typically 
sufficient  for  the  log-normal  approximation  to  be  made. 

PACS  numbers:  43.30.Re,  43.30Vh,  43.30Wi,  43.30.Xm  [MBP] 


INTRODUCTION 

As  long  ago  as  World  War  II,  marine  physicists  under 
the  National  Defense  Research  Committee1  observed  that 
natural  disturbances,  such  as  underwater  turbulence  and 
passing  surface  or  internal  gravity  waves,  often  place  the 
ocean-acoustic  waveguide  in  such  a  state  of  flux  that  a  sig¬ 
nal,  deterministic  when  transmitted  from  the  source,  be¬ 
comes  fully  randomized  after  propagating  only  several  chan¬ 
nel  depths  away  in  range  to  a  receiver.  When  not  properly 
accounted  for,  such  randomization  can  severely  degrade  the 
accuracy  of  an  experimental  measure,  as  well  as  any  param¬ 
eter  estimates  based  upon  this  measure.  However,  with  the 
knowledge  that  the  field  is  fully  randomized,  and  thereby  can 
be  described  by  a  circular  complex  Gaussian  random 
(CCGR)  variable,2  many  useful  statistical  properties  of  sub¬ 
sequent  intensity  measurements  and  parameter  estimates  can 
be  readily  deduced  by  respective  applications  of  coherence 
theory  and  estimation  theory. 

The  CCGR  field  assumption  has  a  long  history.  In  the 
analysis  of  random  signals  and  noise,  it  has  been  legitimately 
made  when  the  central  limit  theorem  applies,  such  as  in  the 
scattering  of  radiation  from  fluctuating  targets3  and  surfaces 
with  wavelength-scale  roughness.4  In  the  present  context,  it 
is  the  basis  for  the  Rayleigh  scatter  channel  that  is  not  only 
frequently  cited  in  communication  theory,5  but  has  also  been 
used  to  describe  the  saturated  region6  of  multipath  propaga¬ 
tion  in  the  ocean  for  many  years.1  Previous  analyses  in  this 
area,  however,  have  been  implicitly  limited  to  certain  special 
cases  for  which  the  time-bandwidth  product  of  the  field  re¬ 
ceived  from  a  given  source  is  unity.  In  this  paper,  the  statis¬ 
tical  formulation  is  extended  and  generalized  to  be  a  function 
of  measurement  time  and  temporal  coherence.  This  more 
general  framework  is  therefore  highly  relevant  to  modern 
ocean-acoustic  sonar  and  communication  systems  that  em¬ 
ploy  time-bandwidth  products  exceeding  unity,  or  average 


over  many  independent  samples.  Such  averaging  is  often  es¬ 
sential  when  the  goal  is  to  reduce  the  variance  of  a  measure¬ 
ment,  or  subsequent  parameter  estimate,  enough  for  the  error 
to  fall  within  a  tolerable  threshold. 

Much  of  the  literature  in  ocean-acoustic  transmission 
scintillation6-11  closely  follows  or  is  based  upon  the  classic 
analysis  of  Dyer,12  who,  like  his  predecessor  Bergmann,1 
concluded  that  in  many  circumstances  a  Rayleigh  scatter 
channel  can  be  used  to  describe  the  random  fluctuations  of 
measured  intensity.  In  order  to  extend  this  analysis  to  the 
case  of  variable  time-bandwidth  product,  some  specific  re¬ 
sults  from  Dyer’s  “Statistics  of  sound  propagation  in  the 
ocean”  12  are  used  as  a  guide  line.  First,  techniques  well 
established  in  coherence  theory  are  used  to  formulate  the 
statistical  properties  of  averaged  intensity  in  terms  of  mea¬ 
surement  time  and  temporal  coherence.  This  is  done  in  Sec. 
I  where  it  is  also  shown  that  the  “short-time  average”  of 
Dyer  corresponds  to  a  measurement  that  has  a  time- 
bandwidth  product  of  unity  and  is  therefore  by  definition 
instantaneous.  In  Sec.  II,  the  distribution  for  the  logarithm  of 
averaged  intensity  is  given  as  a  function  of  measurement 
time  and  temporal  coherence.  As  a  result.  Dyer’s  well  known 
5.6  dB  transmission  loss  (TL)  standard  deviation  and  2.5-dB 
augmentative  bias  in  mean  TL  are  found  to  be  valid  only  for 
the  special  case  of  an  instantaneous  measurement.  The  gen¬ 
eral  TL  standard  deviation  and  mean  have  forms  that,  respec¬ 
tively,  approach  zero  and  —  101og(the  mean-square  trans¬ 
mission),  as  they  must,  in  the  deterministic  limit  of 
increasing  time-bandwidth  product.  Further  analysis  along 
this  line  indicates  that  the  logarithm  of  averaged  intensity  has 
a  distribution  that  rapidly  converges  to  a  Gaussian  for  time- 
bandwidth  products  exceeding  four.  The  significance  of  this 
particular  finding  is  that,  by  invocation  of  the  central  limit 
theorem  for  sufficiently  large  time-bandwidth  products,  in¬ 
tensity  statistics  in  the  saturated  region  can  legitimately  be 
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described  as  log-normal,  which  runs  counter  to  some  previ¬ 
ous  suggestions.12 

Additionally,  certain  special  assumptions  were  made  in 
the  derivation  of  a  well-known  probability  distribution  for 
the  “noise  of  multiple  distant  sources”  12  that  limits  its  gen¬ 
eral  usefulness.  Specifically,  the  intensity  contribution  of 
each  source  was  assumed  to  be  independent,  exponentially 
distributed,  and  comprised  of  a  single  tone  spectrally  disjoint 
from  the  simultaneously  measured  tones  of  the  other  sources. 
These  assumptions  imply  that  the  time-bandwidth  product  of 
the  total  received  field  must  equal  the  number  of  independent 
sources,  and  the  measurement  period  must  be  such  that  the 
spectral  contribution  of  each  tone  is  resolvable  according  to 
the  Rayleigh  criterion. 13  In  Sec.  Ill  of  the  present  paper,  the 
more  general  assumption  is  made  that  the  measured  fields 
from  the  independent  sources  are  independent  CCGR  vari¬ 
ables,  and  no  restriction  is  made  on  their  spectra  except  that, 
for  practical  considerations,  they  have  finite  bandwidth.  The 
intensity  distribution  for  the  “noise  of  multiple  distant 
sources”  derived  under  this  more  general  assumption  is 
given  as  a  function  of  the  measurement  time  and  temporal 
coherence  of  the  total  received  field.  Similar  differences  be¬ 
tween  the  distributions  for  a  “signal  plus  noise”  12  derived 
in  previous  work  and  those  derived  here  are  discussed. 

The  statistics  of  averages  of  independent  intensity 
samples  are  also  investigated.  Such  averages  are  widely  used 
in  a  variety  of  ocean  acoustic  applications.  For  example,  in 
displaying  the  beamformed  output  of  a  hydrophone  array,  it 
is  common  practice  to  reduce  the  variance  by  averaging  the 
uncorrelated  intensities  received  on  adjacent  nonoverlapping 
beams.14  Similarly,  it  is  sometimes  convenient  to  average 
independent  multipath  arrivals  that  are  temporally  disjoint, 
or  to  average  independent  measurements  of  backscatter  to 
reduce  the  variance  in  scattering  strength  or  target  strength 
estimation. 14-15  In  Sec.  IV,  the  probability  distributions  for 
such  amalgamated  intensity  measurements  are  investigated. 
The  probability  distribution  for  the  difference  between  two 
intensity  measurements  is  then  given  in  Sec.  V  to  address  the 
issue  of  monitoring  a  moving  source  or  scatterer. 

A  brief  discussion  of  classical  parameter  resolution 
bounds  and  Fisher  information  is  then  provided  in  Sec.  VI 
for  measurements  obeying  the  various  distributions  pre¬ 
sented.  This  is  used  to  show  that  parameter  resolution  is 
highly  dependent  upon  the  way  that  intensity  measurements 
are  made.  For  example,  the  logarithmic  measures  commonly 
used  in  ocean  acoustics,  such  as  scattering  strength,  target 
strength  and  TL,  must  be  derived  from  a  corrected  version  of 
the  sonar  equation  that  accounts  for  an  inherent  bias  depen¬ 
dent  on  the  time-bandwidth  product  of  the  intensity  average. 
This  bias  attains  its  maximum  magnitude  of  2.5  dB  for  an 
instantaneous  sample  and  only  vanishes  in  the  deterministic 
limit  of  large  time-bandwidth  product.  The  logarithmic  mea¬ 
sures  then  have  mean-square  errors  that  approximate  the 
Cramer-Rao  lower  bound  with  increasing  accuracy  for  in¬ 
creasing  time-bandwidth  product.  Finally,  a  quantitative 
measure  is  given  for  the  amount  of  information  that  can  be 
lost  by  certain  widely  practiced  procedures  for  reducing  a  set 
of  measurements  to  a  single  mean  statistic.  Such  reduction  is 


often  employed  in  ocean-acoustic  processing  but  may  be  det¬ 
rimental  to  subsequent  parameter  estimates. 

I.  INTENSITY  STATISTICS  AS  A  FUNCTION  OF 
MEASUREMENT  TIME  AND  TEMPORAL  COHERENCE 

While  the  expected  intensity  of  a  CCGR  field  may,  for 
all  practical  purposes,  be  temporally  invariant  during  a  given 
set  of  measurements,  the  experimental  measures  of  intensity 
are  themselves  still  subject  to  statistical  fluctuation  regard¬ 
less  of  the  amount  of  finite-time  averaging  employed.  How¬ 
ever,  under  the  more  stringent  assumption  of  stationarity,2 
many  statistical  properties  of  the  these  fluctuations  become 
invariant  and  can  be  readily  expressed  in  terms  of  the  mea¬ 
surement  averaging  time  and  the  coherence  of  the  received 
field.  As  a  first  step  in  this  direction,  the  distribution  for  the 
average  intensity  measured  from  a  CCGR  field,  along  with 
the  first  few  moments,  are  derived  in  terms  of  the  time- 
bandwidth  product  of  the  received  field.  The  derivation,  pre¬ 
sented  in  Sec.  I  A,  mechanically  parallels  that  given  in  the 
optics  literature16  for  the  statistical  properties  of  polarized 
thermal  light.  The  way  that  this  analysis  extends  previously 
derived  ocean-acoustic  intensity  statistics  is  discussed  in  Sec. 
I  B.  A  brief  summary  of  relevant  literature  in  acoustics,  ra¬ 
dar,  and  optics  is  then  given  in  Sec.  I  C. 

A.  The  gamma  distribution 

Let  the  field  measured  at  a  receiver  be  denoted  by 
z(f)exp(— i2irfct).  The  envelope  z(t)  contains  the  stochastic 
properties  of  the  field  modulated  at  constant  carrier  fre¬ 
quency  f  c .  Let  both  the  modulated  and  demodulated  fields 
be  CCGR  variables  such  that  the  real  x(t)  and  imaginary 
y(t)  components  of  the  envelope  z(t)  are  independent 
Gaussian  random  variables  with  zero  mean,  and  the  same 
variance.  Therefore,  at  any  instant  t,  the  probability  that  the 
envelope  will  have  value  z  =  x  +  iy  is 

1  [  x2  +  y2] 

for  —oo  <x,y<°°,  (1) 

where  the  variance  of  z  is  given  by  (x2)  +  (y2)  and 
(x2)  =  (y2).  Rewriting  Eq.  (1)  in  polar  coordinates,  the  uni¬ 
form  phase  distribution  can  be  integrated  out.  The  squared 
magnitude  of  the  field  then  defines  the  instantaneous  inten¬ 
sity  I(t)  =  \z(t)\2,  which  obeys  an  exponential  distribution 


=  0,  elsewhere,  (2) 

with  mean  I  =  (I)  =  2(x2)  and  variance  (/)2  =  (/2)  — (7)2.  An 
essential  quality  of  this  exponential  form  is  that  the  most 
probable  value  for  instantaneous  intensity  is  the  same  as  the 
most  probable  value  for  the  field,  namely  zero. 

While  the  concept  of  an  instantaneous  intensity  will 
presently  be  shown  to  be  more  than  a  mere  mathematical 
stepping  stone,  it  certainly  must  be  used  in  this  way  to  for¬ 
mulate  the  statistical  properties  of  actual  intensity  measure¬ 
ments.  This  is  because  actual  intensity  measurements  can 


770  J.  Acoust.  Soc.  Am.,  Vol.  100,  No.  2,  Pt.  1,  August  1996 


Nicholas  C.  Makris:  Saturated  transmission  scintillation  770 


only  be  made  over  a  finite  time  period  but  still  need  to  be 
expressed  in  terms  of  the  instantaneous  fluctuations  of  an 
underlying  random  field.  For  example,  let  the  average  inten¬ 
sity  measured  over  time  interval  T  be 

1  rr/2 

W=-  I(t)dt.  (3) 

‘  J  -TI2 


A  measure  of  the  number  of  independent  intensity  fluctua¬ 
tions  or  coherence  cells  /i  averaged  during  the  period  T  is 
given  by  the  squared-mean-to-variance  ratio  or  signal-to- 
noise  ratio  (SNR)  of  the  measurement  W.  Under  the  assump¬ 
tion  that  the  field  undergoes  stationary  fluctuations,  such  that 
its  autocorrelation  (z*(t)z(t  +  t))  is  only  a  function  of  the 
time  lag  r  and  not  the  absolute  time  reference  t,  the  SNR  of 
W  is  expressible  in  terms  of  the  temporal  coherence  function 


{z*(t)z(t+  t)> 

7(r)'  <Rop> 


(4) 


and  the  triangle  function 


elsewhere. 


(5) 


which  arises  from  auto-convolution  of  a  rectangular  mea¬ 
surement  window  of  length  T.  Specifically,  the  SNR  of  W,  or 
equivalently  the  number  of  independent  intensity  fluctuations 
averaged  during  the  measurement  time  T,  is  not  restricted  to 
discrete  values  but  rather  is  defined  by  the  continuous  vari¬ 
able 


W2  = 
^  {w2)-{w)2 


1 


(6) 


where  the  last  equality  is  obtained  using  CCGR  moment 
factoring.2 

Apart  from  experimental  intrusion,  the  random  field  un¬ 
dergoes  fluctuations  that  occur  over  a  characteristic  period  tc 
referred  to  as  its  coherence  time  scale.  A  useful  measure  of 
the  coherence  time  scale  is 


(7) 


Noting  the  infinite  integration  limits  in  this  definition,  one 
would  expect  that  more  accurate  experimental  estimates  of  tc 
would  be  obtained  for  longer  measurement  times  T.  This  is 
certainly  the  case,  and  can  be  easily  shown  by  considering 
the  limiting  form  of  Eq.  (6)  for  T>  tc  .  Here,  the  number  of 
coherence  cells  is  well  approximated  by  the  linear  relation¬ 
ship  /i=TI tc  between  the  measurement  period  and  coher¬ 
ence  time  scale.  In  the  opposite  extreme  of  a  measurement 
time  much  less  than  the  field’s  intrinsic  coherence  time  scale, 
such  that  T<tc,  the  integral  of  Eq.  (6)  approaches  its  mini¬ 
mum  value  of  unity,  and  /x=  1  becomes  a  very  good  approxi¬ 
mation.  In  this  case,  the  average  intensity  measurement  W  is 
nearly  instantaneous,  and  for  practical  purposes  obeys  the 
exponential  distribution  of  Eq.  (2),  but  not  in  the  vicinity  of 
zero  unless  [jl  is  identically  equal  to  one,  as  will  be  demon¬ 
strated  presently. 
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FIG.  1.  Time-bandwidth  product  /i  as  a  function  of  T/tc  ,  where  T  is  the 
measurement  period  and  rc  is  the  intrinsic  coherence  time  scale  of  the  fluc¬ 
tuating  field.  The  solid  line  shows  this  dependence  when  the  field  has  a 
rectangular  frequency  spectrum.  The  dashed  line  shows  the  linear  asymptote 
for  large  T/rc . 


First,  it  is  useful  to  consider  the  general  dependence  of  /x 
on  T/tc  as  illustrated  in  a  particular  example.  This  depen¬ 
dence  is  given  in  Fig.  1  for  the  case  in  which  z(t)  has  a 
rectangular  frequency  spectrum  such  that  ■':/(/  )  =  tc  for 
|/|  =£  1/(2  rc)  and  •9e’(f)  =  0  elsewhere,  )  being  the 

Fourier  transform  of  M  r)  which  takes  the  form 


y(  r)  = 


sin(7 tt/tc) 

77  t!  Tc 


(8) 


It  is  easy  to  verify  that  in  an  estimate  of  the  spectrum  of  z(t) 
over  finite  time  window  T,  only  j±  independent  frequency 
components  can  be  resolved  by  the  Rayleigh  criteria  in  the 
interval  |/|=£1/(2tc),  which  equals  the  number  of  temporal 
samples  attainable  at  the  Nyquist  rate.17  In  this  case,  as  in  the 
general  case  of  an  integrable  spectrum,  fi  measures  the  time- 
bandwidth  product  of  the  field  received  over  finite-time  win¬ 
dow  T,  just  as  the  ratio  /jJT  measures  its  bandwidth.  For 
example,  when  /i  is  linearly  dependent  upon  T,  as  it  is  for 
T>  tc  ,  the  measured  bandwidth  fitT  is  a  good  approxima¬ 
tion  to  the  intrinsic  bandwidth  of  the  fluctuating  field  1  Itc  .  In 
the  opposite  extreme  when  tc>T ,  the  measured  bandwidth 
of  roughly  1  IT  is  so  dominated  by  the  effect  of  temporal 
windowing  that  it  cannot  provide  useful  information  about 
the  coherence  time  scale  of  the  fluctuating  field.  (A  more 
rigorous  discussion  of  the  time-bandwidth  product  can  be 
found  in  Ref.  18.) 

In  this  context,  if  the  continuous  time  average  W  is  re¬ 
placed  by  an  ensemble  average  of  /a  independent  and  iden¬ 
tically  distributed  instantaneous  intensity  samples  from  the 
same  data,  both  the  mean  and  variance  of  the  resulting  en¬ 
semble  average  are  identical  to  those  of  the  original  continu¬ 
ous  average.  Furthermore,  the  probability  distribution  for  the 
ensemble  average  is  readily  obtained  as  the  inverse  Fourier 
transformation  of  the  characteristic  function2  for  instanta¬ 
neous  intensity  raised  to  the  power  /i.  This  process  leads  to 
the  well-known  gamma  distribution 
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P(W)  = 


exp {-/x(WII)} 
T(/a) 


for  0^fK<°°, 


=  0,  elsewhere. 


(9) 


that  generally  provides  an  excellent  approximation  to  the  ex¬ 
act  distribution  for  averaged  intensity  W ,  which  has  a  much 
less  convenient  form,  as  is  discussed  in  detail  in  Ref.  2.  For 
example,  Eq.  (9)  converges  to  the  exact  distribution  for  W  in 
both  the  limits  of  T>  tc  and  T<tc.  However,  some  slight 
discrepancy  between  this  approximate  form  and  the  exact 
form  exists  when  T  is  approximately  equal  to  tc  . 

When  instantaneous  intensity  is  digitally  sampled  at  the 
Nyquist  rate,  as  is  common  practice,  Eq.  (9)  is  the  exact 
distribution  for  discrete-time-averaged  intensity.  Equation 
(9),  therefore,  is  essential  to  the  present  analysis  because  it 
provides  a  distribution  for  fluctuating  intensity  in  terms  of 
two  quantities  that  may  be  estimated  with  great  accuracy  and 
relative  ease  in  an  experimental  situation,  namely  the  mean 
intensity  I  and  the  time-bandwidth  product  p.  of  the  measure¬ 
ment.  Additionally,  because  the  average  intensity  W  has  a 
variance  (I)2/p  that  depends  upon  the  mean  /,  its  fluctua¬ 
tions  depend  on  the  expected  value  of  the  signal  and  there¬ 
fore  fall  under  the  category  of  signal-dependent  noise. 

Returning  again  to  the  issue  of  measuring  instantaneous 
intensity  in  an  actual  experiment,  upon  inspection  of  Eq.  (9), 
it  becomes  apparent  that  j±=  I  is  the  only  value  for  p  that 
does  not  exclude  the  possibility  of  measuring  an  average 
intensity  W  that  can  be  identically  zero.  For  all  other  values, 
including  p  very  near  to,  but  not  identically  unity,  the  instan¬ 
taneous  intensity  passes  through  but  does  not  remain  at  the 
isolated  value  of  zero  during  the  measurement  time.  The 
positive-semidefinite  nature  of  instantaneous  intensity  then 
insures  that  a  finite-time  average  can  never  yield  an  average 
value  that  is  identically  zero.  The  forgoing  may  explain  why, 
in  many  experimental  situations  of  high  resolution  and  large 
sample  population,  finely  binned  histograms  of  fluctuating 
intensity  may  follow  the  exponential  distribution  very 
closely  except  in  the  vicinity  of  the  origin,  where  the  fre¬ 
quency  of  samples  decays  to  zero  rather  than  growing  to  a 
maximum  value. 

To  investigate  the  asymptotic  properties  of  the  gamma 
distribution  for  large  time-bandwidth  product,  it  is  conve¬ 
nient  to  employ  the  general  notation 

ffl/(#)  =  <[#-»ii(#)]'),  (10) 

for  the  central  moments  of  integer  order  z>  1 ,  where 
ni  |(i)j=(i)),  and  the  random  variable  ii  obeys  an  arbitrary 
probability  distribution.  It  is  also  convenient  to  employ  the 
compact  notation  ■'£•(/;; ,  ( W),p)  to  denote  the  gamma  distri¬ 
bution,  and  ,/F  {ni\(  rj)jn2(rj))  to  denote  the  Gaussian  distri¬ 
bution 


p  (  1  /  _  [  rj—ni  | ( zy ) ] 2 \ 

v  V  [27tzh2(  ?7)]1/2  6XP\  2 m2(y)  )’ 

where  — oo<^<oo. 

While  it  is  not  difficult,  with  the  aid  of  a  digital  com¬ 
puter,  to  computationally  verify  that  for  large  p.  the  gamma 
distribution  '-0{m ,  ( W),/i)  tends  toward  a  Gaussian 


./V\m  |  (W),m2(W))  with  the  same  mean  and  variance,  an 
analytic  proof  has  been  supplied  by  Mandel,  who  demon¬ 
strated  the  asymptotic  convergence  of  all  the  corresponding 
moments.16  A  compelling  but  less  complete  demonstration 
can  be  made  by  comparison  of  the  kurtosis  (peakedness). 


^  ffl4(#) 

^=[m2(#)]2’ 

and  skew 


(12a) 


nidi}) 

Jd=[m2(tf)]3/2’  (12b) 

of  these  distributions,  where  the  necessary  moments  are 
readily  obtained  upon  differentiation  of  the  appropriate  char¬ 
acteristic  functions.  For  the  Gaussian  distribution,  this  pro¬ 
cedure  leads  to  a  kurtosis  j#v=3,  and  a  skew  jr;=0,  the  latter 
indicating  symmetry  in  the  distribution  of  probability  about 
the  mean.  For  the  gamma  distribution,  the  kurtosis 


d  3  H — , 


and  skew 


(13a) 


2 

4W=~j=,  (13b) 

V/r 

have  asymptotically  Gaussian  behavior  for  p.>\.  or  roughly 
p^  10  in  practice.  However,  averaged  intensity  is  positive 
semidefinite,  and  therefore  its  distribution  never  spans  the 
full  domain  of  a  Gaussian,  except  arguably  in  the  mathemati¬ 
cal  limit  of  arbitrarily  large  pi. 

B.  Application  to  ocean-acoustic  transmission 

As  noted  during  its  development,  the  above  formulation 
for  the  statistical  properties  of  average  intensity  is  based 
upon  the  assumption  that  the  underlying  acoustic  field  at  the 
receiver  undergoes  CCGR  fluctuations.  Due  to  the  recur¬ 
rence  of  CCGR  fields  in  such  a  wide  variety  of  otherwise 
unrelated  realizations  of  stochastic  wave  propagation,  which 
follows  from  the  central  limit  theorem,  it  is  often  difficult  to 
determine  the  cause  of  uncertainty  in  a  particular  saturated 
measurement  without  soliciting  external  information.  In 
ocean  acoustics,  such  external  information  has  been  histori¬ 
cally  collected,  and  a  set  of  basic  causes  for  field  fluctuations 
can  be  readily  cited. 

Specifically,  the  randomization  of  acoustic  fields  in  the 
ocean  typically  arises  from  (1)  incoherent  source  fluctuation, 
(2)  fluctuation  of  a  scatterer,  (3)  relative  motion  between  a 
scatterer  and  source  or  receiver  in  a  waveguide  (4)  relative 
motion  between  source  and  receiver  in  a  waveguide,  (5)  fluc¬ 
tuation  of  the  waveguide  boundary,  as  for  example  due  to 
surface  waves,  (6)  medium  scintillation  or  fluctuation  in  the 
index  of  refraction  due  to  such  phenomenon  as  internal 
waves  or  turbulence.  Typically  the  last  three  are  associated 
with  transmission  scintillation,  which  becomes  saturated  af¬ 
ter  relatively  short  propagation  ranges,  greater  than  a  wave¬ 
guide  depth,  if  many  more  than  a  single  acoustic  mode  con¬ 
tributes  significantly  to  the  received  field,  in  keeping  with  the 
central  limit  theorem.  The  origin  of  the  fluctuation  in  these 
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cases  is  the  change  in  the  interference  structure  of  the  wave¬ 
guide  brought  on  by  the  given  relative  motions,  which  there¬ 
fore  typically  need  to  exceed  only  a  wavelength  scale  in 
random  amplitude  for  saturation  to  occur.  For  example. 
Dyer12  has  noted  that  his  model  of  transmission  scintillation 
“applies  to  an  ensemble  of  experiments  that  contains  signifi¬ 
cant  rangewise  motion  of  the  source  with  respect  to  the  (mul¬ 
tipath  interference)  structure,”  as  in  case  (4)  above.  How¬ 
ever,  it  is  the  point  of  this  paper  that  such  motion  contributes 
to  the  temporal  fluctuations  whose  time  scale  must  be  com¬ 
pared  to  the  measurement  averaging  time  to  obtain  general 
and  accurate  statistics.  Even  early  researchers  in  ocean- 
acoustics  found  that  these  fluctuations  can  occur  over  ex¬ 
tremely  short  periods,  less  than  1  min,1  so  highlighting  the 
need  for  the  more  general  formulation  provided  here. 

The  validity  of  the  CCGR  field  assumption  can  easily  be 
tested  in  practice  by  examining  the  statistics  of  the  field  re¬ 
ceived  by  a  hydrophone  over  time.  The  extension  to  previous 
work  lies  in  the  realization  that  the  temporal  coherence  func¬ 
tion  or  spectrum  at  the  receiver  should  also  be  estimated  to 
determine  the  time-bandwidth  product  of  the  measurement. 
This  then  leads  to  the  general  distribution  for  averaged  in¬ 
tensity  as  a  function  of  measurement  time  and  temporal  co¬ 
herence  given  here  in  Eq.  (9). 

In  this  context,  the  exponential  distribution,  given  in  Eq. 
(9)  of  Ref.  12  for  intensity  measured  in  a  “short-time  aver¬ 
age,”  is  only  valid  for  measurements  of  instantaneous  inten¬ 
sity  where  p,=  1 ,  and  not  for  longer  stationary  averaging  pe¬ 
riods  for  which  /i>  1 .  While  Dyer  defines  a  “short-time 
average”  as  “...  an  average  taken  over  a  time  long  as  com¬ 
pared  to  lirho  (the  carrier  frequency’s  period)  but  much  less 
than  the  stationarity  time  ...”,12  this  definition  is  not  suffi¬ 
cient  to  meet  the  necessary  requirement  for  the  exponential 
intensity  distribution  to  be  valid,  namely  that  fi=  1.  If  the 
“coherence  time”  were  substituted  for  the  “stationarity 
time”  in  this  definition,  it  would  be  stated  correctly  in  the 
present  context. 

The  contention  of  Ref.  12  that  the  variance  of  the  natural 
log  of  averaged  intensity  is  “independent  of  any  of  the  met¬ 
rics  of  the  problem,”  therefore,  must  be  put  into  perspective. 
This  is  because  the  metric  describing  the  number  of  coher¬ 
ence  cells  /jl  in  the  intensity  average  was  implicitly  assigned 
to  unity.  The  relevance  of  these  issues  to  the  well  known  5.6 
dB  transmission  loss  standard  deviation  of  that  reference  will 
be  examined  in  Sec.  II. 

C.  Historical  notes 

Almost  immediately  after  Bergmann’s  analysis  of 
ocean-acoustic  intensity  fluctuations  for  the  war  effort,  and 
much  earlier  than  its  unclassified  appearance  more  than  20 
years  later.  Rice19  showed  that  the  gamma  distribution  could 
be  used  to  describe  the  statistical  properties  exhibited  by  the 
finite  time  average  of  an  exponentially  distributed  random 
variable.  Shortly  thereafter,  Mandel16  applied  Rice’s  work  in 
his  demonstration  that  the  time-integrated  intensity  of  polar¬ 
ized  thermal  light  undergoes  fluctuations  that  can  be  well 
described  by  the  gamma  distribution,  which  has  lain  the 
foundation  for  the  formulation  given  in  Sec.  I  A  of  the 
present  paper. 


With  a  minor  change  of  variables,  the  same  distribution 
was  obtained  independently  by  Swerling  in  his  case  II  prob¬ 
ability  for  the  radar  detection  of  a  fluctuating  target  by  a 
mean-square  processor  which  integrates  over  returns  from  /jl 
active  pulses.3  An  ocean-acoustic  confirmation  of  Swerling’ s 
work  has  apparently  been  given  in  the  underwater  target 
scattering  histograms  of  Dahl  and  Mathisen.20  Those  histo¬ 
grams  were  derived  by  peak  amplitude  analysis,  correspond¬ 
ing  to  a  time-bandwidth  product  of  unity,  for  which  the  am¬ 
plitude  distribution  is  Rayleigh  and  the  intensity  distribution 
is  exponential.  Similar  concepts  and  distributions  have  been 
used  to  describe  the  intensity  statistics  of  laser  beam  fluctua¬ 
tions  induced  by  propagation  through  the  turbulent  atmo¬ 
sphere.  In  that  community  the  model  is  known  as  the  ‘  ‘single 
scatter  model”  of  atmospheric  scintillation.21'22 

Images  derived  from  CCGR  fields  exhibit  the  same  kind 
of  signal-dependent  noise  described  by  the  statistical  formu¬ 
lation  given  in  Sec.  I  A.  The  noise  in  such  images  is  com¬ 
monly  referred  to  as  speckle  when  the  time-bandwidth  prod¬ 
uct  of  the  measurement  is  near  unity,  for  which  the  intensity 
fluctuations  attain  their  largest  variance.  Such  noise  is  often 
found  in  synthetic  aperture  radar  (SAR),  medical  ultra  sound, 
side-scan  sonar  and  towed-array  reverberation  images  be¬ 
cause  the  associated  active  systems  typically  irradiate  terrain 
of  wavelength-scale  roughness  with  narrow-band  waveforms 
of  low  time-bandwidth  product. 

II.  STATISTICS  OF  LOGARITHMIC  INTENSITY 
MEASURES  AND  TRANSMISSION  LOSS 

It  is  traditional  in  ocean-acoustics,  and  many  other  dis¬ 
ciplines,  to  measure  fluctuating  intensity  in  logarithmic  units, 
although  the  reasons  given  for  such  a  measure  often  seem  to 
verge  on  scientific  folklore.  For  example,  reference  is  often 
made23  to  the  apparent  logarithmic  response  of  human  audi¬ 
tory  and  visual  perception  to  intensity  stimulus,24  the  impli¬ 
cation  being  that  such  a  response  has  been  optimized  by 
millions  of  years  of  evolution.  However,  a  particularly  com¬ 
pelling  quantitative  advantage  of  the  logarithmic  measure 
exists  and  can  be  readily  cited.  That  is,  for  the  average  in¬ 
tensity  of  a  fully  randomized  Gaussian  field,  a  logarithmic 
measure  homomorphically  transforms  signal-dependent 
noise  into  additive  signal-independent  noise.  Consequently, 
optimal  methods  for  finding  expected  signals  or  patterns  in 
independent  additive  noise,  which  also  happen  to  be  well 
established,  can  be  directly  applied  when  the  fluctuating  in¬ 
tensity  of  a  CCGR  field  is  measured  in  logarithmic  units. 

In  this  context,  a  theoretical  foundation  for  logarithmic 
measures  of  fluctuating  intensity  has  recently  been  derived  in 
the  optics  literature.25  It  is  based  upon  criteria  for  optimal 
pattern  recognition  that  follow  from  statistical  estimation, 
optimal  filter  and  information  theory.  Because  this  theory 
supports  the  use  of  traditional  logarithmic  measures  of  fluc¬ 
tuating  intensity,  it  should  not  be  surprising  that  the  probabil¬ 
ity  distributions  upon  which  it  is  based  also  describe  the 
statistics  of  ocean-acoustic  TL  measurements.  This  formula¬ 
tion  is  presented  here  in  Sec.  II  A.  The  way  it  extends  pre¬ 
vious  work  is  discussed  in  Sec.  II  B.  A  brief  summary  of 
relevant  literature  in  acoustics  and  optics  is  given  in  Sec. 
II  C. 
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In  a  tangential  point,  the  question  of  whether  the  theory 
of  Ref.  25  has  any  bearing  on  the  logarithmic  response  of 
human  auditory  and  visual  perception  to  intensity  stimulus 
has  been  recently  raised,26  primarily  because  the  optical  and 
acoustics  fields  received  by  the  eye  and  ear  often  undergo 
CCGR  fluctuations. 


A.  The  exponential-gamma  distribution 

The  intensity  level  is  defined  by 


W 

L  =  In—,  (14) 

7  ref 


where  /ref  is  a  reference  intensity  such  that  the  normalized 
mean  intensity  is  70  =  ///ref.  Following  Bayes’  theorem,  the 
gamma  distribution  for  the  measurement  W  leads  to  the 
exponential-gamma  distribution 


P(L)  = 


(/*/ 1  exp  {- p(eL/I0)  +  L} 

T  (p) 


for  — °°<L<  oo,  (15) 

for  the  log-transformed  statistic  L ,  which  is  denoted  by 
<?8?({L),p)  for  notational  convenience.  The  central  mo¬ 
ments  of  the  exponential-gamma  distribution  are  most 
readily  obtained  by  direct  integration  with  respect  to  the  vari¬ 
able  W.  Such  a  process  yields  the  expected  intensity  level 

(L)  =  ln(7o)  +  iA(^)-ln /r,  (16a) 

and  variance 

(L2)-(L)2=a  2,(1),  (16b) 

which  in  general  is  not  inversely  proportional  to  the  time- 
bandwidth  product  p.  Here,  i/Ap)  is  Euler’s  psi  function27 
and  £(2,/z)  is  Riemann’s  zeta  function.27  Formulas  defining 
these  special  functions  are  given  in  Appendix  A.  For  ex¬ 
ample,  £(2,1  )  =  7r/6,  and  in  the  limit  as  p>  1,  the  expectation 
value  of  L  converges  to  ln(/0)  while  the  variance  £(2 ,p)  ap¬ 
proaches  l//x. 

It  is  of  great  practical  significance  that  intensity  level 
L  has  a  variance  that  does  not  depend  on  the  expected  inten¬ 
sity  /,  as  average  intensity  W  does,  but  only  on  the  time- 
bandwidth  product  p  of  the  measurement.  For  example,  sup¬ 
pose  that  the  instantaneous  intensity  of  a  CCGR  field  is 
practically  stationary  over  short  periods,  but  not  over  long 
periods  where  some  trend  in  expected  intensity  emerges. 
Further  suppose  that  this  trend  is  to  be  measured  by  time 
series  analysis,  where  average  intensity  samples  are  consecu¬ 
tively  collected  over  short-term  stationary  periods  and  then 
concatenated.  The  resulting  intensity  time  series  will  have  a 
standard  deviation  that  is  directly  proportional  to  the  local 
value  of  expected  intensity,  making  comparison  of  samples 
with  different  expectation  values  difficult.  Conversion  to 
logarithmic  units,  however,  homomorphic  ally  transforms 
such  signal-dependent  noise  into  additive  signal-independent 
noise.  The  log-transformed  time  series  has  a  uniform  stan¬ 
dard  deviation  when  the  time-bandwidth  product  of  the  mea¬ 
surement  is  chosen  to  be  constant  for  all  samples.  Well- 
established  techniques  for  processing  signals  in  additive 
noise  are  then  appropriate.  In  particular,  the  optimal  method 


FIG.  2.  (a)  Kurtosis  for  the  exponential-gamma  distribution  in  solid  line  and 
for  the  gamma  distribution  in  dashed  line.  The  Gaussian  kurtosis  is  3.  (b) 
Skew  for  the  exponential-gamma  distribution  in  solid  line  and  for  the 
gamma  distribution  in  dashed  line.  The  Gaussian  skew  is  zero. 


for  recognizing  the  nonstationary  trend  is  provided  by 
matched  filtering  the  time  series  with  hypothetical  trends  in 
the  log-transformed  domain.25  This  technique  has  also 
proven  to  be  valuable  for  recognizing  patterns  in  active  sonar 
images.14'28  The  speckle  found  in  such  images,  which  is  most 
pronounced  for  low  time-bandwidth  product  measures,  arises 
from  the  same  CCGR  field  fluctuations  described  here.  Char¬ 
acteristic  scales  of  the  nonstationary  trends  in  both  saturated 
intensity  time  series  and  speckled  images  are  best  measured 
in  the  logarithmic  domain.25 

Some  of  the  asymptotic  properties  of  the  exponential- 
gamma  distribution  can  be  ascertained  by  considering  its 
kurtosis 


and  skew 

~  [£(2,/r)]3/2 ’ 


(17a) 

(17b) 


for  increasing  p.  Apparently,  the  exponential-gamma  ap¬ 
proaches  a  Gaussian  with  the  same  mean  and  variance  more 
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rapidly  than  the  gamma  distribution,  as  is  evident  after  in¬ 
spection  of  the  curves  provided  in  Fig.  2.  For  practical  pur¬ 
poses,  the  criterion  that  /x>4  appears  to  be  sufficient  for 
eFS<((L),/r)  to  be  approximated  by  yV({L),{L2)  —  (L)2). 

The  Gaussian  asymptote  of  the  exponential-gamma  dis¬ 
tribution  can  also  be  obtained  by  analytic  means.  After  some 
straight  forward  algebraic  manipulation,  Eq.  (15)  can  be  re¬ 
written  as 

^(^)  =  (lY— jj  exp(/r{[L-ln(70)]-exp[L-ln(70)]}). 

(18) 

The  argument  of  the  interior  exponential  can  then  be  ex¬ 
panded  into  a  Taylor  series  such  that 

P(L)~(l7/0  exp<-/*))exp(  “f  lL_lnUo)l2J. 

(19) 

for  |L-ln(/0)|^3.25'29  This  can  be  used  as  a  relatively  crude 
approximation  even  for  /n=  1  because  the  corresponding 
standard  deviation  of  L  is  then  7r/v6  =  1.3.  As  //  increases  the 
approximation  becomes  better.  By  use  of  Stirling’s  formula29 
for  /i>  1 ,  Eq.  (19)  transforms  to  the  Gaussian 

P(L)~[j^)  exp(-f  |£-ln(/o)|2),  (20) 

so  that  the  exponential-gamma  distribution  for  the  intensity 
level  L  is  distributed  according  to  .  'V'{\n(I {^),  l//x)  for  large 
//.  Therefore,  as  the  time-bandwidth  product  of  the  measure¬ 
ment  becomes  large,  the  expected  intensity  level  approaches 
the  logarithm  of  expected  intensity,  the  intensity  level  vari¬ 
ance  approaches  the  inverse  time-bandwidth  product,  and  the 
exponential-gamma  distribution  converges  to  a  Gaussian. 

B.  Application  to  ocean-acoustic  transmission  loss 
statistics 

Due  to  its  simple  linear  dependence  on  intensity  level  L, 
transmission  loss  //,  sampled  after  saturated  ocean-acoustic 
propagation,  has  statistical  properties  that  can  be  readily  ex¬ 
pressed  in  terms  of  the  time-bandwidth  product  of  the  mea¬ 
surement.  Specifically,  TL  is  related  to  intensity  level  by  the 
equation 

H=-l0log(eL)  +  K,  (21) 

where  H  is  measured  in  dB  re:  1  m,  and  K  is  a  conversion 
factor  in  dB  re:  1  //Pa  and  m  that  can  be  set  to  zero  without 
any  loss  of  generality.  Therefore,  the  probability  distribution 
for  H  is  readily  found  to  be  exponential-gamma,  so  that  the 
mean  TL 

10  log  e)[ln(7o)  + </<//)  — In  //],  (22) 

not  only  depends  on  the  mean  intensity  /,  but  also  on  the 
time-bandwidth  product  of  the  measurement.  The  TL  stan¬ 
dard  deviation 

(H2)-(H)2=(  10  log  e )sqrt[£(2,//)],  (23) 

on  the  other  hand,  only  depends  upon  the  time-bandwidth 
product.  When  these  moments  are  plotted  as  a  function  of  //, 
as  in  Fig.  3,  it  becomes  evident  that  Dyer’s  5.6  dB  TL  stan- 


FIG.  3.  Transmission  loss  (TL)  standard  deviation  sqrt[(//2)  —  {H)2]  in 
solid  line  and  TL  mean  ( H )  for  /()=1  in  dashed  line  as  a  function  of  time- 
bandwidth  product  fjL.  The  asymptotic  dependence  of  the  TL  standard  de¬ 
viation  for  large  //,  (10  log  e)sqrt[  l///j.  is  also  given  in  dot-dashed  line. 

dard  deviation  and  2.5-dB  augmentative  bias  in  mean  TL  are 
only  valid  for  an  instantaneous  measurement,  for  which 
fi=  1,  and  not  for  longer  stationary  averages  as  might  be 
inferred  from  Ref.  12.  As  the  time-bandwidth  product  in¬ 
creases,  the  TL  standard  deviation  asymptotically  approaches 
zero  along  the  curve  (10  log  e)sqrt[l//z],  or  roughly 
4.34sqrt[l//z],  while  the  TL  mean  approaches  —10  log  (the 
mean-square  transmission)  as  it  must  in  the  deterministic 
limit  of  an  arbitrarily  large  sample  size.  Additionally,  the 
asymptotic  analysis  of  the  previous  section  shows  that  TL 
statistics  in  the  saturated  region  can  legitimately  be  described 
as  normal  for  measurement  time-bandwidth  products  exceed¬ 
ing  four,  which  runs  counter  to  some  previous  suggestions 
implicitly  based  upon  unity  time-bandwidth  measures.12 


C.  Historical  notes 

What  is  referred  to  here  as  the  exponential-gamma  dis¬ 
tribution  was  apparently  first  derived  by  Dyer12  for  an  en¬ 
tirely  different  purpose,  namely  to  describe  the  “noise  of 
multiple  distant  sources.”  Flowever,  the  relationship  be¬ 
tween  degrees  of  freedom,  integration  time  and  temporal  co¬ 
herence  was  not  addressed,  as  is  discussed  further  in  the  next 
section.  The  distribution  was  later  rediscovered  by  Barakat4 
in  an  analysis  of  laser  beam  speckle  patterns  observed 
through  a  finite  spatial  aperture.  In  Barakat’ s  analysis,  the 
relationship  between  degrees  of  freedom,  spatial  aperture  and 
spatial  coherence  is  derived  in  accord  with  Goodman.30  It  is 
noteworthy  that  Pierce  has  arrived  at  a  decibel  standard  de¬ 
viation  equivalent  to  4.34sqrt(l//z)  dB  for  Gaussian  random 
acoustic  signals  without  exploiting  the  convenience  of  com¬ 
plex  variables.31  The  derivation  of  the  Gaussian  asymptote  of 
the  exponential-gamma  distribution  follows  the  analysis  of 
Arsenault  and  April29  and  Makris.25 
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III.  THE  AVERAGE  INTENSITY  OF  MULTIPLE 
INDEPENDENT  CIRCULAR  COMPLEX  GAUSSIAN 
RANDOM  FIELDS 

When  multiple  independent  CCGR  fields  are  superposed 
at  a  receiver,  fluctuations  in  average  intensity  can  once  again 
be  described  by  the  gamma  distribution,  as  is  shown  in  Sec. 
Ill  A.  The  way  that  this  extends  and  generalizes  previous 
statistical  formulations  for  the  “noise  of  multiple  distant 
sources”  and  a  “signal  plus  noise”  is  discussed  in  Sec. 
Ill  B. 

A.  The  gamma  distribution 

When  multiple  independent  fields  are  superposed  at  a 
receiver,  their  individual  variances,  or  equivalently  their  ex¬ 
pected  intensities,  sum  to  the  variance  of  the  total  field, 
which  constitutes  the  expected  intensity  at  the  receiver. 
However,  the  intensities  of  the  independent  fields  may  not 
comprise  mutually  independent  statistics  of  the  measurement 
except  under  certain  special  circumstances. 

To  illustrate  this  situation,  let  the  total  received  field 
z(f)exp(— i2  7 rfct)  be  the  sum  of  5  independent  CCGR  fields 
Zj(f)exp(— Hirffif,  emanating  from  distinct  sources  or  scat¬ 
tered,  modulated  at  carrier  frequency  fc  such  that 
5 

z(0  =  2  Zi(t).  (24) 

i  =  1 

Then  it  is  readily  verified  with  characteristic  functions 
that  z(t)  is  a  CCGR  variable  with  variance  1=  (|z2|)  equal  to 
the  sum  of  the  variances  of  the  zft),  which  are  denoted 
by  Ij  =  (\zf\)-  The  expected  intensity  of  the  superposed  fields 
is  then 

7=2  f.  (25) 

i  =  1 

While  the  average  intensity  W  of  the  superposed  fields,  mea¬ 
sured  during  period  T,  still  obeys  the  gamma  distribution 
•5-(/,/x),  the  time-bandwidth  product  /x,  as  defined  in  Eq.  (6), 
is  now  dependent  upon  the  temporal  coherence  function  of 
the  superposed  fields.  Furthermore,  when  the  variance  of  W , 
namely  (/)2//r,  is  expanded  in  terms  of  the  intensities  of  the 
component  fields,  cross  terms  with  factors  /,/,  emerge.  These 
cross  terms  appear  because  the  instantaneous  intensities  of 
the  constituent  fields  are  not  independently  superposed  in  an 
instantaneous  measurement  of  the  total  field  intensity.  An 
alternative  proof  of  this  is  given  in  Appendix  B. 

Suppose  now  that  the  desired  signal  Zi(t),  with  inten¬ 
sity  /| ,  is  just  one  component  among  many  received  in  a 
noisy  environment,  where  the  noise  field  is  z(f)  —  Z\(t),  with 
intensity  IN=I—Ii.  Given  the  average  intensity  measure¬ 
ment  W,  the  SNR  for  the  intended  signal  component  then 
becomes  /xf/j//)2,  which  includes  the  effects  of  both  addi¬ 
tive  noise  and  signal-dependent  held  fluctuations.  Here,  the 
time-bandwidth  product  fi  depends  upon  the  coherence  func¬ 
tion  of  the  combined  signal  and  noise  fields.  Consequently,  if 
the  additive  noise  extends  over  a  much  broader  frequency 
band  than  the  signal,  stationary  averaging  will  reduce  the 
additive  noise  component  of  the  variance  of  W  far  more 
rapidly  than  the  signal-dependent  component.  A  more  effec¬ 


tive  way  of  eliminating  out-of-band  noise  is  to  filter  the  total 
received  held  to  the  signal  band,  as  is  typically  accomplished 
with  a  matched  hlter. 


B.  Application  to  ocean-acoustic  transmission 
statistics 

The  distributions  analyzed  in  Sec.  Ill  A  of  this  paper, 
namely  Eqs.  (9)  and  (15)  with  I  defined  for  multiple  sources, 
also  characterize  the  “noise  of  multiple  distant  sources” 
studied  by  Dyer. 12  However,  there  are  substantial  differences 
between  Dyer’s  formulation  and  that  given  here.  Specifi¬ 
cally,  the  distributions  defined  in  Eqs.  (27)  and  (28)  of  Ref. 
12  were  derived  under  the  assumption  that  the  intensity  mea¬ 
sured  from  each  distant  source  is  independent,  exponentially 
distributed  and  comprised  of  a  single  tone  spectrally  disjoint 
from  the  simultaneously  measured  tones  of  the  other  sources. 
However,  the  condition  of  Ref.  12  that  each  source  be  of 
“different  frequency”  is  not  sufficient  to  insure  that  the  in¬ 
tensity  contribution  measured  from  any  one  source  will  be 
independent  of  that  measured  from  any  other.  The  additional 
requirement  for  such  independence  is  that  the  inverse  mea¬ 
surement  time  1/T  must  be  less  than  the  spectral  separation 
between  any  of  the  contributing  fields.  Further,  the  intrinsic 
bandwidth  of  the  field  received  from  each  source  must  be  far 
less  than  1/T  for  the  measured  contribution  from  that  source 
to  be  virtually  instantaneous,  as  is  necessary  for  its  intensity 
distribution  to  be  approximately  exponential.  These  require¬ 
ments  can  be  stated  mathematically  by  expressing  the  time- 
bandwidth  product  in  terms  of  the  spectrum  S^(fi  )  of  the 
total  received  field 


sin(TrT(f—f')\2 

TTif-f)  I 


df  df 


1 


(26) 

which  is  simply  another  representation  of  Eq.  (6).  The  first 
integration  is  a  convolution  which  sets  the  number  of  inde¬ 
pendent  frequency  cells,  while  the  second  integration  sums 
these  cells.  The  cells  are  measured  in  units  of  1/T,  the  dis¬ 
tance  from  the  apex  of  the  sine  function  to  its  first  zero 
crossing,  in  accord  with  the  Rayleigh  Criterion.  Therefore, 
the  previous  statistical  description  of  “the  noise  of  multiple 
distant  sources,”  based  upon  the  assumption  that  the  inten¬ 
sities  received  simultaneously  from  the  noise  sources  are 
independent,12  comprises  a  special  case  of  the  more  general 
formulation  provided  here,  where  it  is  assumed  that  the  fields 
received  from  the  noise  sources  are  independent. 

The  distributions  for  “signal  plus  noise”  of  the  present 
section  differ  substantially  from  those  given  in  Eqs.  (33)  and 
(35)  of  Ref.  12  for  reasons  similar  to  those  just  discussed. 
Here,  for  example,  TL  for  a  “signal  plus  noise”  field  obeys 
the  exponential-gamma  distribution  with  standard  deviation 
(10  log  e)sqrt[£(2,/u)]  dB,  where  the  time-bandwidth  product 
/jl  is  for  a  measurement  of  the  combined  fields  with  expected 
intensity  I=IX  +  IN.  Therefore,  the  TL  standard  deviation  for 
a  “signal  plus  noise”  field  never  exceeds  5.6  dB  as  it  may  in 
Ref.  12. 
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IV.  THE  AVERAGE  OF  INDEPENDENT  GAMMA- 
DISTRIBUTED  INTENSITY  MEASUREMENTS 


In  ocean-acoustics,  there  are  many  situations  in  which 
independent  and  presumably  gamma-distributed  intensity 
measurements  are  averaged.  For  example,  in  displaying  the 
beamformed  output  of  a  hydrophone  array,  it  is  common 
practice  to  reduce  the  variance  by  averaging  the  uncorrelated 
intensities  received  on  adjacent  nonoverlapping  beams.14 
Similarly,  it  is  sometimes  convenient  to  average  independent 
multipath  arrivals  that  are  temporally  disjoint,14  or  to  average 
independent  measurements  of  backscatter  to  reduce  the  vari¬ 
ance  in  scattering  strength  or  target  strength  estimation. 15  To 
investigate  the  statistical  properties  of  such  measurements,  a 
distribution  is  derived  for  the  average  of  two  independent 
gamma-distributed  intensity  samples  that  may  not  be  identi¬ 
cally  distributed.  The  asymptotic  Gaussian  form  of  this  dis¬ 
tribution  is  then  generalized  to  describe  the  average  of  an 
arbitrary  number  of  independent  but  not  necessarily  identi¬ 
cally  distributed  intensity  samples.  This  form  is  later  used,  in 
Sec.  VI,  to  quantify  the  maximum  amount  of  information 
that  can  be  inferred  about  a  desired  parameter  set  from  an 
amalgamated  measurement. 

Let  the  average  intensities  Wa  and  Wp  be  independent 
and,  respectively,  obey  the  nonidentical  gamma  distributions 
'^•(/,,,/x,,)  and  '£(1  p, /i p).  The  distribution  for  the  sum 
k=W a+W p  is  readily  obtained  as  the  inverse  Fourier  trans¬ 
form  of  the  product  of  the  characteristic  functions  of  W a  and 
Wp.  With  the  aid  of  tabulated  integral  transforms,27'32  the 
distribution  for  k  is  found  to  be 


(Kfla/Ia)tJ-^KfJ,p/Ip)fl^  [Lp  I 

P(  K)  =  -  expj  —  -=T-  K 

KY(lLa+Hp)  [  Ip  j 


X  \F i 


Up  P*  a 

Ha\Ha+ Hp\K\  -f —  y 
W/3  *a  , 


for  /c^O, 


(27) 


and  zero  elsewhere,  where  1F1  is  Rummer’s  confluent  hyper¬ 
geometric  series.32  The  mean  of  k  is  (k)=I a+ 1 p,  and  the 
variance  is 


<*>-(*>  = 


da)2  ,  dp)2 


P-a 


P*  P 


(28) 


with  a  minor  change  of  variables,  this  equation  describes  the 
probability  distribution  for  partially  coherent  thermal 
light.2'33  The  mean  of  k  is  still  (k)=I a+ 1 p  while  the  vari¬ 
ance  simplifies  to 


( x2)-(x)2  = 


da)2  +  dp)2 

P- 


The  kurtosis 


(30) 


,  da)4+dp)4 

\W  \da)2+dp)2f’ 


and  skew 


(31a) 


_/  2  \  da)3  +  dp )3 

w ^hdS+dp)2?12 


(31b) 


indicate  that  Eq.  (29)  converges  to  the  Gaussian 
lT^/t),^2)— (k)2)  for  1,  as  a  natural  consequence  of  the 
central  limit  theorem. 

When  it  is  the  average  of  the  two  measurements  that  is 
of  interest,  the  random  variable  k  must  be  redefined  as 
K=(Wa+Wp)/ 2.  This  average  is  distributed  according  to 
Eq.  (29)  but  with  Ia  replaced  by  I J2  and  Ip  by  I  J.  Assum¬ 
ing  the  statistics  across  the  two  measurements  are  indepen¬ 
dent  and  ergodic,  such  that  I=Ia=Ip,  the  average  k  obeys 
the  distribution  m  ,2/x),  as  is  consistent  with  a  doubling  of 
the  time-bandwidth  product  of  a  constituent  sample.  By  as¬ 
ymptotic  analysis  of  Eq.  (29),  a  rigorous  demonstration  of 
this,  and  the  case  when  Ip  approaches  zero,  is  provided  in 
Appendix  C.  Similarly,  for  an  ergodic  population  of  indepen¬ 
dent  samples  distributed  according  to  3?(/,/r),  the  average  of 
D  samples  of  this  population  obeys  ■5(1  J)  /a). 

The  more  general  scenario  is  the  ensemble  average  of  D 
independent  intensity  measurements  W,  that,  respectively, 
obey  the  nonidentical  gamma  distributions  SW->, /*,■).  Let 
this  be  described  by 

1  D 

k=j (32) 

U  i=\ 


For  Hj>  1,  the  probability  distribution  for  k  converges  to  the 
Gaussian  (k)2),  where  respective  linear  sum¬ 

mations  of  the  constituent  means  and  variances  lead  to  the 
mean 


It  is  often  possible  to  make  the  simplification  that  the  time- 
bandwidth  product  of  the  two  measurements  is  equal  so  that 
j±(=  p,p= ji.  The  distribution  for  the  sum  k,  readily  obtained 
by  convolution  of  the  distributions  for  Wa  and  Wp,  is  then 

P[K)  =  Jj_  / !±  V~ 
r (fi)\Ta-Tff)  U  TJp  I 


<*>=i2  (wd, 

D  i=  i 


and  variance 

,  !w  ,2  1  y  W2 

(*>-<«> 


(33a) 


(33b) 


flK  da  +  Ip)  [  /  flK_  dg-Ip)  \ 

2  Tjp  j  m_i/2\  2  rjp  r 


for  a:  3®  0,  (29) 

and  zero  elsewhere,  where  I  l/2  is  a  modified  Bessel  func¬ 
tion  of  the  first  kind  of  order  /r— 1/2.  It  is  noteworthy  that. 


of  the  average,  as  follows  from  the  independence  of  the  con¬ 
stituent  measurements.  If  now  each  W t  represents  the  aver¬ 
age  intensity  of  .S’,  independent  CCGR  fields,  the  expectation 
value  for  each  W,  is 

Si 

(W,)  =  E  Iu,  (34) 

7=1 
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in  accord  with  the  analysis  of  the  previous  section.  The  util¬ 
ity  of  this  representation  will  become  evident  in  Sec.  VI  and 
Appendix  E,  where  a  quantitative  measure  of  the  amount  of 
information  that  can  be  lost  by  such  an  averaging  process, 
common  in  ocean-acoustic  applications,  is  presented. 


V.  THE  RATE  OF  CHANGE  OF  AVERAGE  INTENSITY 


It  is  sometimes  possible  to  infer  the  velocity  of  a  moving 
source  or  scatterer  from  a  series  of  incoherent  intensity  mea¬ 
surements  or  images.34  When  this  is  possible,  it  is  necessary 
to  compute  the  time  rate  of  change  of  measured  intensity  by 
the  method  of  finite  differences.  While  such  a  velocity  esti¬ 
mate  may  also  require  computation  of  finite  differences  in 
space  to  account  for  advection  in  the  image  plane,  the  dis¬ 
tinction  between  temporal  and  spatial  domain  is  not  crucial 
to  the  present  statistical  analysis. 

Suppose  two  consecutive  averaged  intensity  measure¬ 
ments  Wa  and  W p  are  separated  by  a  time  interval  At  long 
enough  to  insure  that  they  are  independent.  Let  these  mea¬ 
surements  obey  respective  nonidentical  gamma  distributions 
•'5-(/ri,/x)  and  '-0(1  p,f±).  The  random  variable  for  intensity 
rate  is  then  defined  by  the  difference  quotient 
V=(W p—  Wa)/At.  Its  probability  distribution  is  readily  ob¬ 
tained  by  inverse  Fourier  transformation  of  the  relevant  char¬ 
acteristic  functions.  With  the  aid  of  the  tabulated  integral 
transforms,27  the  distribution  for  V  is  found  to  be 


—  ,  —  ,  \  1/2 


P(V)= 


1  U|V|A tyhiAt  (Ia+Ip) 


F(Al)  \  Ia  +  h 


IJ/3 


Xexpj  — 


fiAt  \V{Ip-Ia)\ 


IJ/3 


fl\V\\t  ( Ia  +  Ip)'' 


IJp 


for  —  00<  y<oo 


(35) 


where  Ky 2-^  is  a  modified  Bessel  function  of  the  second 
kind  of  order  1/2— /jl.  The  mean  is  (V)  =  (I  p— I  a)  /  At  and  the 
variance  is 


(v2)-(v)2= 


2_{Ia)2  +  (Ip)2 


ma  ty 


(36) 


The  kurtosis //y ,  given  by  the  right-hand  side  of  Eq.  (31a), 
and  skew 


_2_ 

y  1 « [(/j2+(/W]”- 

indicate  that  Eq.  (35)  converges  to  the  Gaussian 
V),(V2)  —  (  V)2)  for  /jl>  1,  in  accord  with  the  central 
limit  theorem. 

To  obtain  the  second  time  derivative  of  intensity.  A, 
three  measurements  are  needed.  If  these  are  equally  spaced 
in  time,  the  distribution  of 


A  = 


Wa-2Wp+Wy 
(A  t)2 


(38) 


converges  to  the  Gaussian  yK((A),(A2)  —  (A)2)  for  fi>l, 
with  mean 


<A>  = 


Ia~'2-IP+Iy 

(A  t)2 


and  variance 


(39a) 


(A2)-(A)2  = 


(Ia)2+(2Ip)2+(Iy)2 

MA  f)4 


(39b) 


To  compute  higher-order  derivatives,  well-known  finite  dif¬ 
ference  equations  can  be  used.35  However,  the  variance  is 
proportional  to  the  weighted  sum  of  the  variances  of  the 
constituent  measures,  and  gains  positive  terms  linearly  with 
the  order  of  the  derivative  to  be  estimated.  Therefore,  for 
fixed  /jl ,  the  SNR  of  the  derivative  estimate  typically  de¬ 
creases  as  the  order  of  the  derivative  increases. 


VI.  THE  RESOLUTION  OF  PARAMETERS 


Parameter  resolution  is  highly  dependent  upon  the  way 
that  intensity  measurements  are  made.  To  show  this,  classi¬ 
cal  parameter  resolution  bounds  and  Fisher  information  ma¬ 
trices  are  derived  for  the  various  kinds  of  measurements  ana¬ 
lyzed  in  previous  sections.  In  particular,  the  logarithmic 
measures  commonly  used  in  ocean  acoustics,  such  as  scat¬ 
tering  strength,  target  strength  and  TL,  must  be  derived  from 
a  corrected  version  of  the  sonar  equation  that  accounts  for  an 
inherent  bias  dependent  on  the  time-bandwidth  product  of 
the  intensity  average.  The  logarithmic  measures  then  have 
mean-square  errors  that  approximate  the  Cramer-Rao  lower 
bound  (CRLB)  with  increasing  accuracy  for  increasing  time- 
bandwidth  product.  Additionally,  it  is  shown  how  informa¬ 
tion  can  be  lost  by  certain  widely  practiced  procedures  for 
reducing  a  set  of  independent  samples  to  a  single  mean  sta¬ 
tistic,  as  is  often  done  in  ocean-acoustic  processing. 

Suppose  that  a  general  Na  -dimensional  parameter  vector 
a  is  to  be  estimated  from  the  A'-dimensional  measurement 
vector  Y.  According  to  estimation  theory,36  the  mean-square 
error  of  any  unbiased  estimate  a,,  based  upon  measurement 
vector  Y,  can  never  be  less  than  the  CRLB36 

E[(at-ai)2]Hr1(*)h,  (40) 


where  at  is  the  true  parameter  value  and  J(a)  is  the  Na  by  Na 
Fisher  information  matrix36  with  elements 

[  d2  1 


Jlj(a)  =  -E 


da  j  daj 


In  P(Y|a) 


(41) 


A.  Information  in  intensity  and  log-transformed 
intensity  measures 


Suppose  the  parameters  a  are  to  be  estimated  from  a  set 
of  N  independent  gamma-distributed  intensity  measurements 
W) r.  contained  in  the  vector  W.  Because  no  information  is  lost 
in  the  homomorphic  transformations  Lk=ln(Wk/ITef),  where 
the  Lk  comprise  the  vector  L,  the  Fisher  information  matri¬ 
ces  for  measurements  W  and  L  are  identical  and  equal  to25 

Y  P'k  ^t(®)  dl] fc(®)  ,  , 

JiJ(a)=2J  -= - -r - ,  (42) 
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as  shown  in  Appendix  D.  This  expression  exhibits  behavior 
that  should  be  expected  in  a  measure  of  information. 
Namely,  information  is  positive  semidefinite,  and  the  infor¬ 
mation  contained  in  a  joint  set  of  independent  measurements 
equals  the  cumulative  information  of  all  the  individual  mea¬ 
surements.  Consequently,  the  CRLB  on  parameter  resolution 
is  only  singular  when  the  Fisher  information  for  all  joint 
measurements  is  zero,  as  is  evident  by  inspection  of  the  par¬ 
ticular  case  given  in  Eq.  (42). 

Although  both  intensity  measurements  and  their  log- 
transforms  contain  the  same  information,  the  information 
about  parameters  linearly  related  to  their  respective  expecta¬ 
tion  values  is  substantially  different.  For  example,  given  the 
single  gamma-distributed  measurement  W,  the  CRLB  for  es¬ 
timation  of  I  is  simply  the  variance  of  the  measurement 
J~ 1  (/)  =  (I)1!  /jl,  so  that  the  measurement  W  itself  attains  the 
CRLB  as  an  unbiased  estimate  of  I. 

The  situation,  however,  is  different  for  logarithmic  mea¬ 
sures.  By  straightforward  manipulation  of  Eq.  (42),  the 
Fisher  information  for  the  measurements  W  or  L  can  be 
rewritten  as25 


Jij(  a) 


d(Lk)  d{Lk) 
2j  AL  ~ L - t: 


(43) 


with  the  understanding  that  the  fik  do  not  depend  on  a.  Given 
the  single  exponential -gamma-distributed  measurement  L, 
the  CRLB  for  estimation  of  ( L )  therefore  is  J  ]  ( (L) )  =  \  / /i. 
While  the  variance  of  L  is  approximately  equal  to  this  bound 
even  for  fi=l,  it  only  converges  to  it  in  the  asymptotic  limit 
of  large  /jl ,  as  shown  in  Sec.  II.  The  CRLB,  therefore,  can 
only  be  attained  asymptotically  for  parameters  linearly  re¬ 
lated  to  the  expected  value  of  intensity  level.  One  should 
realize,  however,  that  for  practical  purposes  such  asymptotic 
convergence  occurs  for  relatively  small  /jl,  as  noted  in  Sec. 
II. 


B.  The  corrected  sonar  equation 

If  the  estimate,  based  upon  the  exponential-gamma- 
distributed  measurement  L,  is  for  the  unknown  decibel  quan¬ 
tity 

a  =  (10  log  e)[ln  7q  +  fl0],  (44) 

as  it  often  is  in  practice,  where  fl0  is  some  known  constant, 
the  unbiased  estimator 


The  familiar  sonar  equation  a  =  (10  log  e)[L+fl0], 
therefore,  is  only  valid  in  the  deterministic  limit  of  large  /jl. 
The  corrected  sonar  equation  given  in  Eq.  (45)  includes  the 
additional  terms  necessary  to  account  for  the  /r-dependent 
bias  introduced  by  an  exponential-gamma-distributed 
intensity-level  measurement  and  is  therefore  the  more  gen¬ 
eral  form. 


C.  Information  in  the  average  intensity  of  superposed 
fields 


The  situation  becomes  more  complex  when  each  Wk 
measures  the  average  intensity  of  a  superposition  of  Sk  inde¬ 
pendent  CCGR  fields,  as  in  Sec.  Ill  where  reception  of  the 
“noise  of  multiple  distant  sources”  and  a  “signal  plus 
noise”  is  considered.  In  this  case,  the  Fisher  information 
matrix  of  Eq.  (42)  becomes 


A/O) 


N 


2 


zSk 

r—  1 


dlqk{&)  dlrk( a) 
da,  daj 


[tfiM  a)]2 


(46) 


This  form  is  useful  because  it  explicitly  shows  how  informa¬ 
tion  partitioned  among  the  component  fields  is  incorporated 
in  an  intensity  measurement  of  their  combination.  Addition¬ 
ally,  it  should  be  understood  that  in  such  a  measurement, 
some  information  carried  by  each  individual  field  is  lost. 
This  is  shown  explicitly  in  Appendix  E. 


D.  Information  in  amalgamated  intensity  averages 

Fisher  information  for  measurements  obeying  more 
complicated  distributions,  such  as  those  described  in  Secs. 
IV  and  V,  is  most  readily  obtained  by  the  methodology  of 
Ref.  37.  Specifically,  suppose  the  parameters  a  only  depend 
on  the  expectation  values  M=(Y).  Then,  given  expressions 
for  M  and  J(M),  the  Fisher  information  matrix  J(a)  can  be 
obtained  from  the  equation38 

c)Mt  o>M 

(47) 

Along  these  lines,  consider  the  amalgamated  average  k  of  D 
independent  gamma-distributed  measurements  W l  that  each, 
respectively,  sample  the  intensity  of  S t  independent  super¬ 
posed  CCGR  fields,  as  described  in  Sec.  IV.  For  the  joint 
measurements  of  N  such  independent  samples  Kk ,  the  Fisher 
information  matrix  is  given  by 


a  =  (10  log  e)[L  +  Fl0—  i/r(/r)  +  ln  /jl\,  (45) 

only  converges  to  (10  log  e)[L+fI0]  in  the  asymptotic  limit 
of  large  /jl,  where  it  also  attains  the  CRLB  on  root-mean- 
square  error  J~  U2(a)  =(10  log  e)sqrt[l//r].  The  more  fre¬ 
quently  encountered  estimator  (10  log  e)[L+ft0],  however, 
has  a  negative  bias  that  decays  from  —2.5  dB  at  /i  =  1  to  zero 
along  the  curve  (10  log  e)[i/»(yu)— In  /r]  as  the  time- 
bandwidth  product  /jl  increases.  Misguided  use  of  this  biased 
estimator  in  reverberation  analysis,  for  example,  can  lead  one 
to  draw  the  erroneous  conclusion  that  scattering  strength  in¬ 
creases  with  the  time-bandwidth  product  of  the  measure¬ 
ment. 


Dk  Smk  qrk(,&i)  dlmnk(&) 


<k  [S?«/„t(a )]2 

A1 Ik 

(48) 

While  this  expression  is  certainly  cumbersome,  it  is  not  with¬ 
out  practical  value.  In  particular,  Eq.  (48)  not  only  quantifies 
the  information  contained  in  the  kind  of  amalgamated  inten¬ 
sity  measurements  that  are  widely  used  in  ocean  acoustics, 
but  also  expresses  this  information  in  terms  of  constants  of 
both  the  measurement  process  and  constituent  fields.  The 


V  W  Y  Sqk  Y 
N  Zjq=\  Zjr=\Zj 


A;/a)=  2 

k=  1 
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amount  of  information  lost  in  such  an  amalgamation,  how¬ 
ever,  can  be  significant,  as  is  shown  in  Appendix  E. 

For  illustrative  purposes,  consider  the  case  when  only  a 
single  parameter  is  to  be  estimated  from  a  single  measure¬ 
ment,  and  that  parameter  is  either  the  expectation  value  of 
the  measurement  or  is  linearly  related  to  it.  Given  the  amal¬ 
gamated  measurement  k ,  the  parameter  Iu,  representing  the 
expected  intensity  of  a  single  field  component  of  the  respec¬ 
tive  intensity  sample  W , ,  cannot  be  resolved  better  than 

,  -  "  [2flA(a)]2 

J~\I n)  =  2  ~~ - •  (49) 


/=  1 


/*/ 


This  bound  equals  Z)2./~  '((/<•)),  and  therefore  is  proportional 
to,  but  potentially  much  greater  than  the  variance  of  the 
amalgamated  measurement  k.  If  such  amalgamation  cannot 
be  avoided,  accurate  estimation  of  the  intensity  Iu  requires 
the  time-bandwidth  products  of  the  constituent  gamma- 
distributed  intensity  samples  to  be  sufficiently  large  that  In 
greatly  exceeds  sqrt[/~  V/n)] . 

In  many  practical  scenarios  in  ocean-acoustics,  however, 
such  amalgamation  cannot  be  avoided.  Consider,  for  ex¬ 
ample,  a  problem  encountered  in  remotely  imaging  the  ocean 
basin  with  an  active  towed-array  system.14  Returns  from 
shadow-zone  sites,  which  lie  between  convergence  zones, 
generally  arrive  in  time  so  that  no  particular  path  to  the  sea¬ 
floor  makes  a  dominant  contribution  during  a  given  measure¬ 
ment  period.  These  returns  may  also  arrive  simultaneously 
and  consecutively  during  the  measurement,  leading  to  an 
amalgamated  average  of  the  intensities  of  superposed  and 
temporally  disjoint  fields.  In  such  cases,  it  is  often  difficult  to 
obtain  sufficiently  large  time-bandwidth  products  to  resolve 
the  mean  contribution  of  a  particular  seafloor  patch.  Returns 
from  shadow  zone  ranges  can  then  be  adequately  described 
as  clutter  because  they  hinder  presently  known  means  of 
inferring  geomorphological  features  of  the  ocean  basin. 


E.  Information  in  intensity-rate  measures 

Finally,  consider  the  intensity  rate  measurement 
V —  (W p —  Wa)/At,  of  Sec.  V,  where  both  W a  and  Wp  are 
gamma-distributed.  By  application  of  Eq.  (47),  the  Fisher 
information  matrix  for  N  independent  measurements  of  in¬ 
tensity  rate  Vk  is 

o 'Tf)k( a)  dTak( a ) \  /  dTPk( a)  dTa]{ a) \ 
ddj  dcij  j  \  daj  da  j  ) 
[^(a)]2+[7^(a)]2 

(50) 

Evidentially,  the  information  about  some  parameter  a,  con¬ 
tained  in  an  intensity-rate  measurement  Vk  may  vanish  when 
the  expectation  value  of  V k  is  zero  even  if  the  parameter 
could  be  uniquely  determined  from  the  constituent  measure¬ 
ments  Wa  or  Wp  .  A  substantial  amount  of  information, 
therefore,  can  be  lost  when  independent  measurements  are 
reduced  by  a  finite  difference  process.  It  is  easy  to  verify, 
however,  that  no  information  about  the  expected  value  of  Vk 
is  ever  lost  when  a  joint  measure  of  Wa  and  Wp  is  replaced 
by  Vk . 


i V  f^k 

■/;/a)  =  2  — 

k=\ 


VII.  CONCLUSIONS 

Coherence  theory  is  used  to  analyze  the  statistical  prop¬ 
erties  of  ocean-acoustic  intensity  fluctuations  measured  after 
saturated  multipath  propagation.  Previous  analyses  in  this 
area  have  been  implicitly  limited  to  certain  special  cases  for 
which  the  time-bandwidth  product  of  the  field  measured 
from  a  given  source  is  unity.  In  this  paper,  the  intensity 
statistics  of  the  saturated  region  are  extended  and  generalized 
to  be  a  function  of  the  measurement  time  and  temporal  co¬ 
herence  of  the  received  field.  The  resulting  intensity  distri¬ 
butions  are  therefore  highly  relevant  to  modern  ocean- 
acoustic  sonar  and  communication  systems  which  employ 
time-bandwidth  products  that  often  exceed  unity,  or  average 
over  many  independent  samples. 

More  general  expressions  are  also  offered  for  the  “noise 
of  multiple  distant  sources.”  Specifically,  it  is  shown  that  a 
previous  and  well-known  assumption  used  to  describe  this 
noise,  that  the  intensities  received  simultaneously  from  the 
noise  sources  are  mutually  independent,  is  a  special  case  of 
the  more  general  assumption  adopted  in  this  paper  that  the 
fields  received  from  the  noise  sources  are  independent. 

The  statistics  of  averages  of  independent  intensity 
samples  are  then  examined  because  such  amalgamation  is 
commonly  employed  in  a  variety  of  ocean-acoustic  measur¬ 
ing  systems.  The  statistics  of  intensity  rate  measures  ob¬ 
tained  by  finite  difference  are  also  examined  to  address  the 
issue  of  monitoring  a  moving  source  or  scatterer. 

A  brief  discussion  of  classical  parameter  resolution 
bounds  and  Fisher  information  for  the  various  distributions 
encountered  is  provided.  This  is  used  to  show  that  parameter 
resolution  is  highly  dependent  upon  the  way  that  intensity 
measurements  are  made.  For  example,  the  logarithmic  mea¬ 
sures  commonly  used  in  ocean  acoustics,  such  as  scattering 
strength,  target  strength  and  TL,  must  be  derived  from  a 
corrected  version  of  the  sonar  equation  that  accounts  for  an 
inherent  bias  dependent  on  the  time-bandwidth  product  of 
the  intensity  average.  This  bias  attains  its  maximum  magni¬ 
tude  of  2.5  dB  for  an  instantaneous  sample  and  only  vanishes 
in  the  deterministic  limit  of  large  time-bandwidth  product. 
The  logarithmic  measures  then  have  mean-square  errors  that 
approximate  the  Cramer-Rao  lower  bound  with  increasing 
accuracy  for  increasing  time-bandwidth  product. 

Finally,  a  quantitative  measure  is  given  for  the  amount 
of  information  that  can  be  lost  by  certain  widely  practiced 
procedures  for  reducing  a  set  of  measurements  to  a  single 
mean  statistic.  Such  reduction  is  often  employed  in  ocean- 
acoustic  processing  but  may  be  detrimental  to  subsequent 
parameter  estimates. 


APPENDIX  A:  EULER’S  PSI  FUNCTION  AND 
RIEMANN’S  ZETA  FUNCTION 


Euler’s  psi  function  is  defined  by 
1  1 


.27 


i/<yU)=-C-2 


k= o  \  n  +  k  k+  1  /  ’ 


1 


1 


for  real  //, 


=  — C+2  7>  f°r  integer  /r>  1, 

k=  i  k 


=  —  C,  for  /i  =  1, 


(Al) 

(A2) 
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where  C  is  Euler’s  constant, 


C  =  lim  X  7- In  fi  =0.577215...  .  (A3) 

/*-=[*=!  k  ) 

Riemann’s  zeta  function  is  defined  by27 

00  1 

£(gaO=X 

k= 0  (/“  +  «) 

for  v>l,  /rAO,—  1,-2,— 3,...  .  (A4) 

APPENDIX  B:  VARIANCE  OF  THE  INSTANTANEOUS 
INTENSITY  OF  A  SUPERPOSITION  OF 
INDEPENDENT  CCGR  FIELDS 

Following  Sec.  Ill  A,  the  instantaneous  intensity  for  a 
sum  of  independent  CCGR  fields  can  be  written  as 

s  s 

/W=X  z,-(*)X  z/w,  (bi) 

i=l  j=  1 

with  expectation  value 


r(/x+  7) r ( /x )  \  4/q,/^ 

v  f  Ua  +  V' 

Xexpl  — /ta: — = — 


After  applying  Gauss’  multiplication  formula  for  gamma 
functions, 

r(nx)  =  (2n)m^-n)n^-mU  r(*+-),  (C2) 

k= 0  \  n! 

for  the  case  n=  2,  and  transforming  the  sum  k=W a+  Wp  to 
the  average  k=(W a+W /3)/2,  it  is  readily  verified  that  Eq. 
(Cl)  reduces  to 


P(k)  = 


(2/jL/I)2tlK2fl-1  exp{-2  fiK/I} 

r(2/t) 


</)  =  £  X  z,-(*)X  z?(f)  =  X  £[|Zi(0|2],  (B2)  Bessel  function  takes  the  forrn/^) 


L«'=i  ;'=i  J  >'=1 

where  the  last  equality  follows  from  the  independence  of  the 
fields.  The  squared  intensity  has  expectation  value 


when  I=Ia=I g.  It  is  noteworthy  that  Eq.  (C3)  also  corre¬ 
sponds  to  the  probability  distribution  for  unpolarized  thermal 
light.2 

Similarly,  as  77  becomes  sufficiently  large,  the  modified 
Bessel  function  takes  the  form  /„(  rj)  e  n!  y  2  tt  77.  Applying 
this  to  Eq.  (29)  for  the  case  I  I  p  the  asymptotic  distribu¬ 
tion  is  simply  8^(1  These  limits  may  also  be  verified  by 
use  of  characteristic  functions. 


(i2)=e  X  z,wX  z*(oX  z*(o2  z*(i) 

i=  1  7=1  7  A:=l  /=  1 

■  5  5  1  r  5  5 

=  £  X  Z;(*)X  Z*(f)  £  X  z*(0X  z*(f) 


*=1  7=1 


*=1  /=  1 


+£  X  z*(f)X  Z*(G  £  X  z,(i)X  zf(f)  , 


*•■1  j  1 


=  1  7=1 


where  the  second  equality  follows  from  CCGR  moment 
factoring.2  Using  Eq.  (B2),  this  can  be  written  as 


is  \  2 

(72)  =  2  X  £[|z,.(t)|2]J  , 


so  that  the  variance  is  the  square  of  the  sum  of  the  expected 
intensities  of  the  superposed  fields 


is  \  2 

</2>-</)2=^X  E[\Zi(t)\2]\  . 


APPENDIX  C:  LIMITING  FORMS  FOR  THE 
PROBABILITY  DISTRIBUTION  OF  THE  SUM  OF  TWO 
INDEPENDENT  GAMMA-DISTRIBUTED  INTENSITY 
MEASUREMENTS 

When  the  expectation  values  of  the  two  independent  in¬ 
tensity  samples  of  Sec.  IV  are  nearly  equal  Ia^Ip,  the  as¬ 
ymptotic  form  I  v(r/)^(r//2)v/T(v+ 1)  for  small  arguments  77 
may  be  substituted  for  the  modified  Bessel  function  in  Eq. 
(29),  assuming  v  is  not  a  negative  integer.  The  resulting  dis¬ 
tribution  for  the  sum  k=W a+W p  is  then 


APPENDIX  D:  FISHER  INFORMATION  FOR 
INDEPENDENT  GAMMA  AND  EXPONENTIAL-GAMMA- 
DISTRIBUTED  MEASUREMENTS 

As  noted  in  Sec.  VI,  the  Fisher  information  for  the  joint 
gamma-distributed  measurements  W,  obeying  the  condi¬ 
tional  probability  distribution 


£(W|a)=II 

k=l 


wkk  *  exP{^/^BV4(a)} 

Hag) 


is  the  same  as  that  for  the  exponential-gamma  distributed 
measurements  L,  obeying  the  conditional  distribution 


£(L|a)=n 

k=  1 


I  exp(Ljfc)  , 

1  AG  —  +A LkEk 

{  Wa) 

Hag) 


It  is  possible  to  demonstrate  this  by  inspection  of  Eq.  (47) 
with  the  additional  knowledge  that  the  matrices  J((W))  and 
J((L))  have  respective  inverses  that  take  the  form  of  the  as¬ 
ymptotic  covariances  of  W  and  L  in  the  Gaussian  limit  of 
large  /ik .  For  example,  J((W))  and  J((L))  are  A  by  A  diag¬ 
onal  matrices  with  respective  elements 

JiMW))  =  ni/(Ti)2Sij,  and  Jij((L))  =  niSij. 

The  more  straight  forward  approach  prescribed  by  Eq. 
(41)  is  instead  presented.  For  the  Fisher  information  matrix 
given  the  measurements  W,  substitution  of  the  right-hand 
side  of  Eq.  (Dl)  in  to  Eq.  (41)  yields 
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Application  of  the  derivative  operator  with  respect  to  cij 
leads  to 
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dlk 
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Zj 
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dai  - 

Taking  the  derivative  with  respect  to  a,  leads  to 

.  ,  ,  J  ^  nk\lwk\dTk  dik 


'  N 

■  2 

(jS-i) 

k=  1 

(4)  1 

\  4  / 

da,  ddj 


+  (h~wk) 
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Finally,  application  of  the  expectation  value  operator  yields 
Eq.  (42)  which  can  be  rewritten  as 


4y(a)=2  AH 

k=  1 


d  In  4(a)  d  In  4(a) 
da :  da ; 


Here,  Ik( a)  can  be  replaced  by  70,  (a)  =/^.(a)//ref  to  reference 
the  logarithmic  measure  to  physical  units  without  altering  the 
Fisher  information. 

For  the  Fisher  information  matrix  given  the  measure¬ 
ments  L,  substitution  of  the  right-hand  side  of  Eq.  (D2)  in  to 
Eq.  (41)  yields 


d  ^  d  -  exp  (Lk) 

J,M)=~E  —  2  -AH  —  In  70  +  -y - 

daik=  i  da -A  k  7n 


Because  VVt/7ref  equals  exp{  Lk)  by  definition,  and  7ref  does 
not  depend  on  a,  Eq.  (D7)  is  identical  to  Eq.  (D3).  Therefore, 
Eq.  (D7)  reduces  to 


■4y(a)=2  AH 

k=  1 


d  In  70t(a)  d  In  70^(a) 


which  equals  Eq.  (D6)  as  expected. 

It  is  noteworthy  that  Eq.  (D6)  can  be  obtained  directly 
from  the  Gaussian  form  of  7>(W|a)  or  P(L|a)  in  the  asymp¬ 
totic  limit  of  /j.k9>  1,  as  is  consistent  with  the  fact  that  both 
the  intensity  measurement  W  and  intensity  level  measure¬ 
ment  L  contain  the  same  information. 


APPENDIX  E:  THE  LOSS  OF  INFORMATION  IN 
COMBINING  NONIDENTICALLY  DISTRIBUTED  FIELDS 
OR  INTENSITIES 

It  is  well  established  that  the  Fisher  information  matrix 
for  a  joint  set  of  independent  measurements  equals  the  sum 
of  the  Fisher  information  matrices  for  each  individual 
measurement.  ’  For  example,  given  scalar  parameter  a, 
the  conditional  probability  distribution  for  the  joint  measure¬ 
ments  of  5  independent  and  instantaneous  CCGR  fields  z,  is 
given  by  the  product 
s 

n  PZi^i\a).  (El) 


The  corresponding  Fisher  information  is  then 

or(£)  ■  <E2) 

Following  Sec.  VI  C,  the  Fisher  information  for  an  instanta¬ 
neous  measurement  of  the  sum  z  of  the  fields  z,  is 

•  (E3) 

where  this  information  is  unchanged  if  the  sum  z  is  replaced 
by  z/jS.  As  a  consequence  of  the  positive  semidefiniteness 
of  expected  intensity  and  the  presence  of  squared  terms,  the 
following  relationship  must  hold  true 


4  f4_ 

4  da 

Ij  da  J 

i=  1  7=1 


[k  =  l  \Ik  da  I  \l=l  m=  1 


=  0,  (E4) 


where  8 )  •  is  the  Kronecker  delta.  By  straight  forward  alge¬ 
braic  manipulation  of  relation  (E4)  it  can  be  shown  that 

'/jomjA  i  7  sum(  a ) .  (E5) 

Therefore,  the  joint  measurements  of  the  S  fields  z,  contain 
more  information  than  a  single  measurement  of  the  sum  z  of 
these  fields.  When  these  fields  are  identically  distributed,  the 
joint  measurements  contain  precisely  S  times  the  information 
of  the  sum.  More  generally,  dividing  the  left-hand  side  of 
relation  (E4)  by 


2  4  , 

k=  1 


gives  the  information  about  parameter  a  that  is  lost  by  re¬ 
taining  the  sum  rather  than  the  joint  measurements. 

Similarly,  given  scalar  parameter  a,  the  conditional 
probability  distribution  for  the  joint  measurement  of  D  inde¬ 
pendent  gamma-distributed  intensity  measurements  W,  is 
given  by  the  product 

D 

II  PwiW;\a).  (E7) 

i  =  i  ‘ 

The  corresponding  Fisher  information  is  then 


^j°int (")  (W,)2 


w\2 

da 


when  the  time-bandwidth  product  of  each  IV,  equals  j±.  Fol¬ 
lowing  Sec.  VI  D,  the  Fisher  information  for  the  average  k  of 
these  D  gamma-distributed  intensities  W,  is 


^average(a)  zf=1(Wj)2  \  da  S  ^  J  ’ 

The  following  relationship  must  hold  true: 

(y  y  n  -  x  J<wj>*<w‘> 

l.-u-i  ij)\<wi)  da  (wj )  da  I 
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because  the  left-hand  side  of  the  equation  must  be  positive 
semidefinite  due  to  the  squared  terms.  By  straight  forward 
algebraic  manipulation  of  relation  (E10)  it  can  be  shown  that 


/joint!  a)^  ^average!  ^  )  > 


(Ell) 


where  equality  holds  when  the  W,  are  identically  distributed. 
The  joint  measurements  of  the  D  intensities  W, ,  therefore, 
contain  more  information  than  a  single  measurement  of  av¬ 
erage  intensity  k  unless  the  intensities  are  identically  distrib¬ 
uted,  in  which  case  the  joint  measurements  and  average  have 
equal  information.  More  generally,  the  left-hand  side  of  re¬ 
lation  (E10)  divided  by 


4  w 

k^\ i  h 


(El  2) 


gives  the  information  about  parameter  a  that  is  lost  by  keep¬ 
ing  the  average  rather  than  the  joint  measurements. 

Finally,  by  letting  (Wi)  =  Ii ,  D  =  S,  /jl=  1,  and  then  com¬ 
paring  Eq.  (E9)  with  Eq.  (E3),  one  can  see  that  information 
is  always  gained  by  measuring  the  instantaneous  intensities 
of  the  S  independent  fields  separately  and  then  averaging 
them,  rather  than  by  first  measuring  the  intensity  of  the  su¬ 
perposed  fields  and  then  dividing  by  S.  Although  this  gain  of 
information  can  be  substantial,  it  requires  the  ability  to  mea¬ 
sure  the  intensity  of  each  independent  field  Z[  separately, 
which  is  not  always  possible  for  the  reasons  given  in  Sec. 
ill  B. 
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